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^ ABSTRACT 

QQ Low-metallicity galaxies exhibit different properties of the interstellar medium (ISM) com- 

pared to nearby spiral galaxies. Obtaining a resolved inventory of the various gas and dust 
components of massive star forming regions and diffuse ISM is necessary to understand how 
those differences are driven. We present a study of the infrared/submillimeter (submm) emis- 
^sf^ sion of the massive star forming complex N158-N159-N160 located in the Large Magellanic 

^-H Cloud. Combining observations from the Spitzer Space Telescope (3.6-70 yum), the Herschel 

Space Observatory (100-500 jim) and LABOCA (on APEX, 870 jim) allows us to work at the 
best angular resolution available now for an extragalactic source (a few parsecs for the LMC). 
We observe a remarkably good correlation between the Herschel SPIRE and LABOCA emis- 
;J] sion and resolve the low surface brightnesses emission. We use the Spitzer and Herschel data 

to perform a resolved Spectral Energy Distribution (SED) modelling of the complex. Using 
modified blackbodies, we derive an average "effective" emissivity index of the cold dust com- 
ponent jSc of 1.47 across the complex. If /3c is fixed to 1.5 , we find an average temperature 
of ~27K (maximum of ~32K in N160). We also apply the GalHano eTaTli^l ) SED mod- 



X 



elling technique (using amorphous carbon to model carbon dust) to derive maps of the star 
formation rate, the grain temperature, the mean starlight intensity, the fraction of Polycyclic 
Aromatic Hydrocarbons (PAH) or the dust mass surface density of the region. We observe 
that the PAH fraction strongly decreases in the H ii regions we study. This decrease coincides 
with peaks in the mean radiation field intensity map. The dust surface densities follow the far- 
infrared distribution, with a total dust mass of 2.1 x 10^ M© (2.8 times less than if carbon dust 
was modelled by standard graphite grains) in the resolved elements we model. We also find a 
non-negligible amount of dust in the region called "N159 South", a molecular cloud that does 
not show massive star formation. We also investigate the drivers of the Herschel/FACS and 
SPIRE submm colours and find that the submm ratios correlate strongly with the radiation 
field intensity and with the near and mid-IR surface brightnesses equally well. Comparing our 
dust map to Hi and CO observations in N159, we then investigate variations in the gas-to- 
dust mass ratio (G/D) and the CO-to-H2 conversion factor Xm- A mean value of G/D~356 
is derived when using 

Xco = 7x10^0 H2 cm-2 (K km s-^y^ (|Fuk ui et al .|2"009|. If a constant 
G/D across N159 is assumed, we derive a Xco conversion factor of 5.4x10^^ H2 cm"^ (K km 

s~^'>~^ We finallv mnrlel inHiviHiial repinns to analvse variations in the SF.D shane aernss the 
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1 INTRODUCTION 

As potential templates of primordial environments, low-metallicity 
galaxies are keystones to understand how galaxies evolve through 
cosmic time. Studying the Interstellar Medium (ISM) of low- 
metallicity galaxies is a necessary step to get a handle on the in- 
terplay between star formation and the ISM under conditions char- 
acteristic of the early universe. Low-metallicity galaxies also have 
quite different infrared (IR) Spectral Energy Distributions (SEDs) 
than solar or metal-rich environments. For instance, their aromatic 
features are diminished compared to dustier galaxies (Madden et al.| 
|2006[ |Engelbracht et al.||2008| ). The paucity of aromatic features 
is usually attributed to the hardness of the radiation field in low- 
metallicity environments, destructive processes such as supernova- 
driven shocks eff'ects ( [Galliano et al. 2003 2005 ) or delayed injec- 
tion of carbon dust by asymptotic giant branch (AGB) stars ( |Dwek| 
|1998[|GaUiano et al.|2008| ) in dwarf galaxies. More relevant to the 
present study, the SEDs of low-metallicity galaxies often exhibit a 
flattening of their submillimeter (submm) slope or a submm excess 
([BottneretaL 2003 , GaUiano et al. 2003 2005 , M arleau et al.|2006[ 
[Bendo et ar][2006| , Galametz et al.,,2009^ ,2011 ), namely a higher 
emission beyond 500 jum than that extrapolated from IR observa- 
tions and standard dust properties (Milky Way dust for instance). 
The origin of this excess is still highly debated. These results high- 
light the importance of a complete coverage of the thermal dust 
emission of low-metallicity objects to get a handle on the overall 
dust population distribution and properties in these environments. 
Combining gas and dust tracers will also allow us to understand 
how the matter cycles and the star formation processes evolve with 
galaxy properties. 

The Large Magellanic Cloud (LMC) is our nearest low- 
metallicity neighbour galaxy (~ 50kpc; |Feast|1999| ), well studied 
at all wavelengths. This proximity enables us to study in detail the 
physical processes at work in the ISM of the galaxy and to individu- 
ally resolve its bright star-forming regions. These structures can un- 
biasedly be isolated due to the almost face-on orientation (23-37^ 
[Subramanian & Subramaniam|2012] ). Furthermore, the low inter- 
stellar extinction along the line of sight facilitates the interpretation 
of the physical conditions of these star-forming regions compared 
to Galactic studies strongly afl'ected by this extinction. Thus, the 
irregular morphology and low-metallicity of the LMC (~ 1/2 Z©; 
|Pagel,2003j ) make it a perfect laboratory to study the evolution in 
the physical properties of the ISM of galaxies and the influence of 
metal enrichment on their star-forming activity. 

Before the Spitzer Space Telescope (Spitzer) observations, the 
infrared studies on the LMC sufl'ered from a lack of wavelength 
coverage or spatial resolution to quantify crucial physical param- 
eters of the ISM properties such as the equilibrium temperature 
of the big grains, the spatial distribution of the dust grain popu- 
lations or the interstellar radiation field. A study of the extended in- 
frared emission in the ISM of the LMC was performed by [Bernard] 
[et al.[ ( [2008| ) as part of the Spitzer Legacy Program SAGE (Sur- 
veying the Agents of a Galaxy's Evolution, Meixner et al.[[2006| 
project using Spitzer observations. They found disparities between 
the overall shape of the SED of the LMC and that of the Milky 
Way (MW), namely a diff'erent mid-infrared (MIR) shape. They 
also found departures from the linear correlation between the FIR 
optical depth and the gas column density. Using Spitzer FIR spectra 
(A = 52-93yL/m), the studies of [van Loon et al.[ ( f2010b| ) also allowed 
comparisons of compact sources in the LMC and its neighbor, the 
Small Magellanic Cloud (SMC), that has an even lower metallicity. 
Their results indicate that while the dust mass diff'ers in proportion 



to metallicity, the oxygen mass seems to diff'er less. The photo- 
electric eff'ect is indistinguishably efficient to heat the gas in both 
clouds. The SMC finally presents evidence of reduced shielding 
and reduced cooling. 

At submm wavelengths, the LMC was observed by [Aguirre[ 
[et al.[ (^2003) with the TopHat instrument, a balloon-borne tele- 
scope ( Silverberg et al. 2003) from 470 yum to 1.2 mm. They con- 
strain the FIR regime with DIRBE (Diff'use Infrared Background 
Experiment) observations at 100, 140 and 240 yt/m and estimated 
an average dust temperature for the LMC of T=25.0+1.8K. Using 
DIRB E and ISSA (IRAS Sky Survey Adas) observations, |Sakon[ 
[et al. |(^2006) found that the submm emission power law index (of- 
ten referred to as /3) is smaller in the LMC than in the MW and that 
the 140 and 240 jim fluxes seem to deviate from the model predic- 
tions, in particular on the periphery of supergiant shells. This excess 
was modelled by a very cold dust component with temperatures < 
9K, even if their lack of submm constraints prevented them to un- 
biasedly conclude that cold dust was the explanation for the excess. 
Using revised DIRBE, WMAP and COBE maps, [Israel etaLl ( |20T0l ) 
construct the global SED of the LMC to confirm a pronounced 
excess emission at millimeter and submm wavelengths, i.e. more 
emission than expected from a submm slope with /5=2. Diff'erent 
hypotheses for this global excess are tested in |Bot et aL] ( |2010a| , 
[Bot et al.| ( [2010b| ), |Galliano et al.| ( |2011[ ) (both rule out for instance 
the "very cold dust" hypothesis) or Planck Collaboration! ( [201 1| ) 
(that test the "spinning dust" hypothesis). 

We more specifically focus the present study on the N158- 
N159-N160 complex, a group of Hii regions (~ 400pc from North 
to South) located in the LMC, ~ 500pc south of the 30 Doradus 
(30Dor) massive star-forming region (catalogued by ' Henize[1956| ). 
The complex is intensively observed in IR, Hi or CO. We have ob- 
tained observations of the complex with the LAB OCA instrument 
at 870 yum, probing the coldest phases of dust at a resolution of 
19''. We note that this analysis presents the first study of LABOCA 
observation in the LMC. We refer to |Bot et al.| (|2010a|) for a de- 
tailed analysis of LABOCA observations of giant molecular clouds 
in the south-west region of the SMC. In addition to LABOCA data, 
the Herschel Space Observatory (Herschel) now allows us to probe 
the thermal dust emission with a coverage of the SED from 70 yt/m 
up to 500 yum, thus to sample the IR peak and the submm slope 
of nearby galaxies. The whole LMC has been mapped with Her- 
schel as part of the HERITAGE project (HERschel Inventory of 
The Agents of Galaxy Evolution; [Meixner et al.|20T0l Meixner et 
al., submitted to AJ). The good resolution of the instruments on- 
board Herschel favours the exploration of the local properties of 
the ISM with a resolution of Herschel/SFIRE at 500 yt/m similar to 
that of Spitzer/MIFS 160 yt/m (FWHM of the PSF=36''), leading, 
for the LMC, to ISM resolution elements of ~9 pc at SPIRE 500 
yL/m. 

This paper thus presents a study of the IR to submm ther- 
mal emission of star-forming regions in various evolutionary states 
at the best resolution available now. This allows us to investigate 
how properties and physical conditions observed on large scales 
are physically linked to specific environments on a more resolved 
scale. In Section 2, we describe the N158-N159-N160 complex and 
the observations and data reduction of the LABOCA, Spitzer and 
Herschel data. We present the SED modelling techniques we apply 
to the Spitzer and Herschel data on resolved scales in Section 3. 
In Section 4, we analyse the results of our modelling, namely the 
distribution of the dust temperatures, the mean stellar radiation field 
intensities, the fraction of Polycyclic aromatic hydrocarbons (PAH) 
or the star formation rates. We also study the dependence of submm 
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Figure 1. The N158-N160-N159 complex is located in the LMC (~ 500pc south of 30Dor). We show here three-color compositions of the two regions N158 
and N160 produced with images taken with HST. For N158 (left), we use a combination of the narrow-band filters F502N/[OIII] (blue), F656N/HQf (green) and 
F673N/[SII] (red). For N160 (right), we use a composition of the narrow-band filters F502N/[OIII] (blue), F656N/HQf (green) and F487N/I^ (red). LABOCA 
870 jim contours are overlaid on both images. Contours are given for 0.018, 0.06, 0.09, 0.15 and 0.3 Jy beam"! on the left image, 0.15 and 0.3 Jy beam on 
the right image. North is up, east is left. 



colours on various parameters and compare our dust mass map with 
the available gas tracers for the N159 region. A SED modelling of 
individual regions across the complex including the 870 fim data in 
the fitting procedure is described in Section 5. We summarize the 
main conclusions of the analysis in Section 6. 



2 A MULTI- WAVELENGTH DATASET 
2.1 The complex 

The N158-N159-N160 complex (~400pc from North to South) is 
located in the LMC at ~500pc south of the 30Dor complex, the 
most massive star-forming region of the Local Group ( [Kennicutt &| 
|Hodge|1986| ). Figure [T] shows three-color compositions of the two 
regions N158 and N160 produced with images taken by the Wide 
Field Planetary Camera 2 (WFPC2) on board the Hubble Space 
Telescope (HST). Data are retrieved from the Hubble Legacy 
Archive (http://hla. stsci.edu/| Program IDs: 11807 and 8247 for 
N158 and N160 respectively). We overlay LABOCA 870 fim 
contours on both images for comparison between the spatial reso- 
lutions. The HST/WFPC2 observations of N158, N159 and N160 
are discussed in details in,Fleener et al.| ( 2010 1, |Heydari-Malayeri| 
|et al.| ( |1999| ) and |Heydari-Malayeri et al.| ( 2002 1 respectively. 

N158 - This Hii region (Fig.[T] left) possesses an elongated 
structure in the N-S direction. Two OB associations were first 
discovered in the region by |Lucke & Hodge| ( |1970| ) who catalogued 
stellar associations across the whole LMC. One is associated with 
the southern part of N 158, one with the north superbubble. Several 
studies have shown that this superbubble is dominated by a young 
stellar population of 2 to 6 Myr ([Testor & Niemela||1998[ |Dunne| 
let al.|2001| ). 



N159 - Observations of N159 in the K band by |Gatley et al.| 
( |1981| ) led to the discovery of the first extragalactic protostar. 
The stellar content of sub- structures of N159 was analysed in 
more detail by Deharveng et al.|(T992| ) at optical wavelengths and 
[Meynadier et al.| ( 2004| l using JHK photometry. N159 presents 
characteristics of on-going star formation activity such as pro- 
tostars, masers, young massive stars as well as High-Excitation 
Blobs (HEB), namely nebulae with turbulent media and strong 
stellar winds interacting with the ambient ionized gas. The Kuiper 
Airborne Observatory (KAO) has been used to probe the complex 
at far-infrared (FIR) wavelengths for the first ti me (|Werner| 
|eraLl|T978l [Jones et al.|[T986l ). Later, [Jones et al.| ( |2005| ) used 
Spitzer/IRAC observations to follow the protostars in N159 and 
constrain the nature of the upper end of the initial mass function 
(IMF) in the LMC. 

N160 - As in N159, massive star formation is well evolved 
in N160 (Fig.[T] right), with parent clouds that have mostly been 
dissipated. The region is associated with H ii regions, young stellar 
clusters and water masers fLazendic et al.[ [20021 [Oliveira et al.[ 
2006 ). Using JHK photometry, Nakaji ma et al.[ ( [2005| ) showed 
that of the 1 1 optical stellar clusters and associations detected in 
the N159/N160 complex (ages listed by [Bica et ar][1996| ), one 
has an age superior to 10 Myr and belongs to N160, the rest 
are younger than 10 Myr, with 3 clusters showing populations 
younger than 3 Myr. We refer to [Carlson et al. (2012) for a study 
of the physical properties and evolutionary stages of YSOs in 
N160 and to [van Loon et al.[ ( [20T0al ) for a detailed study of FIR 
fine-structure cooling lines of compact sources in the 3 Hii regions. 

The N158-N159-N160 complex belongs to the largest molec- 
ular cloud complex of the LMC called "the molecular ridge" and 
accounts for ~30% of the total molecular mass of the galaxy 
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Figure 2. Images of the N158-N159-N160 star-forming complex observed from 8 to 870 iim. North is up, east is left and coordinates (RA, Dec) are given in 
equinox J2000. We provide the intensity scale in MJy sr~^ for each panel. 
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( [Mizuno et al.|200l] ). Several analyses have suggested a dissipated 
sequential cluster formation within N158, N160 and N 159 as well 
as a large-scale sequential cluster formation over the entire com- 
plex (Bolatto et al. 2000 , Nakajima et aL||2005| ) and more glob- 
ally from the 30Dor complex to the southern region of N 159 where 
the large southward CO ridge of 30Dor begins fC ohen et al.|1988[ 
IFukui et al.1|1999| ). Near-infrared (NIR) to FIR observations with 
2MASS and Spitzer highlighted the paucity of star-formation across 
the ridge, with a star-formation activity distributed in lower lumi- 
nosity regions than in the 30Dor complex ( [Indebetouw et al.|2008| ). 
Several studies have suggested alternative scenarios to the standard 
self-propagating star formation to explain this trend, for instance 
that the sequential star formation could be induced in bow- shocks 
formed at the leading edge of the LMC ( |de Boer et al.|1998| ). 

The N159 Hii region harbours two giant molecular clouds 
(GMCs), known as "N1 59 East" (N 159E) and "N159 West" 
(N159W) labeled in Fig. [2|[ChenetaD( 2009) suggest that the cur- 
rent star formation in the N159E is probably triggered by Hii re- 
gions in expansion in the molecular cloud, while the massive stars 
of N159W are forming spontaneously. Another GMC was observed 
in the south of N159 (N159S). Performing a comparative study be- 
tween [CII], CO, and FIR observations, Israel et al. ( 1996 ) find that 
the clouds N160 and N159S are two opposite extremes with respect 
to N 159 (East and West), with molecular gas in N 160 more photo- 
dissociated than in N159. N159S shows a much lower level of star 
formation than N159W and N159E (B olatto et al.|2000[[Rant^^ 
et al. 2005 ), with only a few very faint diff'use Hii regions and no 
OB associations ( |Chen et al. [20101 . The detection of candidates for 
Herbig Ae/Be stars suggest that the cluster formation could have 
just begun in N159S ( Nakaji ma et al.|2005] ). The atomic and molec- 
ular gas distributions across the complex are discussed in more de- 
tails in Section 4.5. 



2.2 LABOCA observations and data reduction 

The full region of N158, N159 and N160 was mapped with 
the LABOCA instrument operating at 870 fim and installed 
on the Cassegrain cabin of the Atacama Pathfinder Experiment 
(APEX) telescope in the Chajnantor site in Chile's Atacama desert. 
LABOCA is a bolometer array of 295 channels with a field of view 
of the full array of 1 1:4 . The full width half maximum (FWHM) 
of its point spread function (PSF) is 19'.'2 + 0.3. Detectors of 
LABOCA are positioned with a double beam channel separation 
(36?) and the field of view is thus under-sampled. The mapping was 
carried out with a raster of spiral pattern to obtain an homogenous 
fully sampled map. A total of 35.5 hrs of observations were ob- 
tained in August 2008 (Program ID: O-081.F-9329A-2008 - Princi- 
pal Investigator: Sacha Hony). Calibration was performed through 
the observations of Mars and secondary calibrators: PMNJ0403- 
8100, PKS0537-441, B13134, V883-ORI, N2071IR and VY-CMa. 
The atmospheric attenuation was determined via skydips every 
hou£] 

The data reduction was performed with the BoA reduction 
package (BOlometer Array Analysis Software]^ The main steps of 
the reduction of the time-ordered data stream of each channel and 

^ See 'http://www.apex-telescope.org/bolometer/laboca/calibration/ for de- 
tails on the calibration and archive tables of sky opacities and calibration 
factors. 

^ BoA was developed at MPIfR (Max-Planck- Institut fr Radio Astronomy, 
Bonn, Germany), AlfA (Argelander-Institut fiir Astronomic, Bonn, Ger- 



scan are: flat fielding, calibration, opacity correction, dead or noisy 
channels removal, removal of the correlated noise on the global ar- 
ray as well as correlated noise induced by the coupled electronics 
of the detectors (amplifier boxes or cables), flagging of stationary 
points or data taken outside reasonable telescope scanning velocity 
and acceleration limits, ^He temperature drift correction, median 
baseline removal and despiking. Individual reduced scans are then 
co-added. 

The presence of very bright structures often leads to negative 
artefacts on their surrounding. These structures are "created" nu- 
merically during the correlation noise or median baseline steps. It 
is possible to reduce these artefacts using an iterative process dur- 
ing the data reduction. Once data are calibrated, we use our reduced 
image to create a "source model" by isolating the pixels superior 
to a given signal-to-noise (5 for the first iterations and decreasing 
progressively to 2). The model grows automatically to avoid iso- 
lated pixels in the masks. We subtract the source model from the 
data and rerun the reduction pipeline. We add the model at the end 
of the new reduction to obtain a new map. This map is used to 
build a new model, input for the following iteration. The process 
is repeated until the reduction converges. Those iterations lead to 
a significant retrieval of faint extended emission around the bright 
structures. Weight maps are derived during the final mapping step 
from which we derive rms and signal-to-noise maps. The final 45 
35' LABOCA image has a pixel size of 9'.'!, with a final rms of 7.8 
mJy beam"^ 



2.3 Herschel Data 

We obtained Herschel data for the N158-N159-N160 complex from 
the HERITAGE project, a programme dedicated to the observa- 
tions of the two Magellanic Clouds and a part of the Magellanic 
Bridge. We refer to |Meixner et al.| ([20T0| and the HERITAGE 
overview of Meixner et al. (submitted) for a detailed description 
of the observing strategy of the HERITAGE project and the data 
reduction of the data. We provide a short summary here. The LMC 
was mapped in two bands with PACS (Photodetector Array Camera 
and Spectrometer; |Poglitsch et al.|2010| ) at 100 and 160 jjm, with 
respective FWHM of the Point-Spread Functions (PSF) of ~7'.'7 
and ~12''. Data are processed in HIPE 7.0 {Herschel Interactive 
Processing Environment) from Level to Level 1 following the 
standard pipeline described in the PACS data reduction guide. Sev- 
eral steps of baseline subtraction (to take 1/^// noise and drifts with 
time into account) and deglitching (to correct for jumps caused by 
cosmic ray hits) are then applied to the data. PACS 100 and 160 
lim timelines are finally mapped using the PhotProject HIPE proce- 
dure. The LMC was also observed with HerschellSVl^^ (Spectral 
and Photometric Imaging Receiver; Griffin et al. 2010) at 250, 350 
and 500 yt/m, with respective FWHMs of their PSFs of 18'^ 25'' and 
36''. Data were reduced in HIPE 7.0, including additional routines 
to subtract the background, adjust the astrometry, apply deglitch- 
ing steps, remove residuals from temperature drifts or mask dis- 
crepant data. We expect residual foreground emission as pointed 
out by Bernard et al., (2008). Using the HI data cube of the LMC 
restricted to velocity ranges matching those of the Galaxy, [Galliano | 
|et al.| ( [20Tl] ) have quantified this contamination to contribute to ~1 
% of the IR power, i.e. smaller than Herschel flux uncertainties. We 



many), AIRUB (Astronomisches Institut der Ruhr-Universitat, Bochum, 
Germany), and IAS (Institut d' Astrophysique Spatiale, Orsay, France) 
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thus consider the foreground contamination to be neghgible in this 
study. 



2.4 Spitzer Data 

The Spitzer observations were performed as part of the SAGE 
project (Surveying the Agents of a Galaxy's Evolution) using both 
instruments IRAC (InfraRed Array Camera; Fazi o et al.|2004| ) and 



MIPS (Multiband Imaging Photometer; Rieke et aL||2004l. The 



IRAC bands cover 3.6, 4.5, 5.8, and 8 jim with FWHMs of the PSFs 
<T'. The MIPS bands cover 24, 70 and 160 jum with a FWHM of 
the PSFs of 18'' and 40'' respectively. A description of the ob- 
serving strategy and data reduction can be found in |Meixner et al.| 
( [2006} . We also refer to [Jones et al.| (12005^ for a very detailed 
overview of the IRAC observations of N 159 and its various com- 
ponents. Because of their lower resolution (40") compared to the 
Herschel/FACS 160 yum maps (12"), we do not use the MIPS 160 
yL/m maps in this study. Meixner et al. (submitted) indicate a good 
agreement of the MIPS and PACS calibrations in the linear range 
of MIPS. 



2.5 Variation of the thermal dust emission with wavelength 

We present a multi- wavelength composition of the N158-N159- 
N160 region in Fig.[2] with, from top to bottom IRAC 8 jim, MIPS 
24 and 70 yt/m, PACS 100 and 160 yt/m, the three SPIRE bands (250, 
350 and 500 yt/m) and the LABOCA map at 870 yt/m. This figure 
shows how the thermal dust emission is evolving with wavelength. 
The 8 yL/m observation mostly traces the aromatic feature emission 
from PAHs. The MIPS 24 yt/m emission is more compact than that 
of the other bands. Tracing the hot dust, the 24 yum emission is 
indeed confined to the compact star-forming sites. The cool dust 
(10-40K) is now traced by the Herschel PACS and SPIRE observa- 
tions. The LABOCA observations are finally probing for the first 
time the coldest phases of dust at a resolution of 19" (4.6pc for the 
LMC). 

We clearly resolve the emission of the bright star-forming re- 
gions of the complex with LABOCA at 870 yt/m, but also detect and 
resolve the diff'use emission at low surface brightnesses. The distri- 
bution of the 870 jim emission is spatially coherent with that of the 
SPIRE maps, both in the H ii regions and in the diff'use medium. To 
allow spatial comparison between the 870 fim maps and the IRAC 
8 jim map, predominantly tracing PAHs, we convolve the IRAC 8 
yL/m map to the LABOCA resolution using the low-resolution kernel 
created by G. Aniano from the PSFs of IRAC and LABOC^and 
observe a spatial coincidence between the 8yum map and the 870yL/m 
contours (Fig. [3]) already observed by [Haas et al.| ( 2002) for a wide 
range of galaxy types. On one hand, this correlation suggests that 
the PAH excitation could be caused by a widespread distribution of 
stellar FUV-emitters and not only linked with star formation ( |Wha-| 
|ley et al.|2009 ). On the other hand, we expect the 870 fim emission, 
namely the cold dust emission, to predominantly originate from 
largely shielded interiors of molecular clouds. Their surfaces are 
known to produce the PAH features ( [Sauvage et al.|2005] ), which 
could explain their association with the 870yL/m emission. 

We finally note the peculiar behaviour of the southern region 
of N159 (N159S in Fig.|2|. As observed in Fig. |2] and Fig. [3] the 
peak of emission moves from 8 to 870 yt/m, with an off'set of 2.7' 



^ http://www.astro.princeton.edu/~ganiano/Kernels.html for kernels and 
details on their construction. 
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Figure 3. Details on the IR/submm emission in the regions. LABOCA 870 
yum contours are overlaid on the IRAC 8 jim observation with, from top 
row to bottom row N158, N160, N159 and N159S. The IRAC 8 yum map is 
convolved to the LABOCA resolution. Contours are given for 0.018, 0.06, 
0.09, 0.15 and 0.3 Jy/beam. North is up, east is left. 
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toward the west for the 870 jdm emission. We will see further in 
this paper that the peak of the 870 jjm emission (which also domi- 
nates the emission at SPIRE wavelengths) spatially coincides with 
the molecular gas reservoir usually referred to as "N159 South" 
(N159S). This region was indeed observed in CO ( [Bolatto et al.| 
[200T||Miz uno et al.l2010l) and is the brightest [Ci] source of N159 
( [Bolatto et al. 2000 ). We provide details on previous observations 
of gas tracers (Hi, CO etc.) in the complex in Section 4.5. 

2.6 Convolution and background subtraction 

We convolve the IRAC, MIPS, PACS and SPIRE 250 and 350 
jim maps to the lowest resolution available, namely that of SPIRE 
500 fim (FWHM: 36''). We used the convolution kernels devel- 
oped by [Aniano et al.| ( |2011| ), which are rotationally symmetric 
and optimized to avoid high-frequency numerical noise (induced by 
the filtering steps when one works in the frequency domain). The 
LABOCA 870 jim map is convolved to the SPIRE 500 jim resolu- 
tion using the low resolution kernel created by G. Aniano from the 
PSFs of both SPIRE 500 yt/m and LABOCA (see Section 2.5). For 
each image, we estimate the background by masking the emission 
linked with the complex and fitting the distribution of the remain- 
ing pixels by a Gaussian. The peak value of this Gaussian is used 
as a background value per pixel and subtracted from the maps. In 
the following study, we restrict ourselves to resolved elements with 
a signal-to-noise superior to 1-cr in the SPIRE bands. 



3 A RESOLVED SED MODELLING 

We use data from 3.6 to 500 jjm to obtain maps of the dust prop- 
erties within the star-forming complex. Because the data reduc- 
tion of the LABOCA data might still lead to removal of some of 
the faintest diff'use emission, we do not use the 870 jim map as a 
constraint for the resolved modelling. We derive average cold tem- 
peratures on a resolved basis using a two-temperature (two modi- 
fied blackbodies) fitting procedure. We then apply a more realistic 
model ( [Galliano et al.|20lT| in order to investigate the distribution 
of the star formation rate, radiation field properties or dust proper- 
ties (mass, PAH fraction) across the star-forming complex. We de- 
scribe the two resolved SED fitting procedures in this section and 
analyse the maps of the physical parameters we obtained in Sec- 
tion 4. In Section 5, we perform a SED modelling of 24 individual 
regions across the complex (same model) to derive integrated prop- 
erties on those regions and investigate the 870 //m thermal emission 
separately. 

3.1 Modified blackbody fitting 

If we assume that dust grains are in thermal equilibrium with the 
radiation field, we can obtain a simple estimate of the cold dust 
temperature across the complex by fitting a modified blackbody 
(MBB) to the dust thermal emission. It is important to keep in 
mind that in reality, each resolution element contains dust grains 
with varied properties along the line of sight. The temperatures de- 
rived using this method should thus be considered as luminosity- 
weighted average values. We choose to fit our 24 to 500yum data 
with a two-temperature (warm+cold) model. This two-MBB fitting 
is a minimum requirement to account for the contribution of warm 
dust to the MIR emission (and estimate the contribution of warm 
dust and single-photon-heating) to the 70 yt/m emission) and avoid 
an overestimation of the cold dust temperatures. 



We fit our data with a model of the form: 

LM, T^, T,,p,) = A, AT^BM. TJ + A, A'^^BM, T,) (1) 

with Tw and Tc the temperature of the warm and cold component, 
fic the emissivity index of the cold component and By the Planck 
function. The A^^ and A^ are scaling coefficients that account for 
the dust masses of each component. The emissivity index of the 
warm dust component is fixed to 2, a standard approximation 
of the opacity in the |Li & Draine| pOOlf dust models. We con- 
volve the model with the instrumental spectral responses of the 
diff'erent cameras to derive the expected photometry. The fit is per- 
formed using the IDL function MPCURVEFIT ( |Markwardt|2009l 
Levenberg-Marquardt least- squares fit). We take the uncertainties 
on flux measurements into account to weight the data during the 
fitting procedure (normal 1/error^ weighting). 

Choice of emissivity - In order to derive an estimate of the 
emissivity index value of the complex, we first let the emissivity 
index of the cold dust component vary {pc free) and perform 
the two-MBB fitting for each ISM resolved elements. We obtain 
an average emissivity index of 1.47 + 0.17 across the complex, 
consistent with the results derived for the whole LMC by Planck] 
[Collaboration! ( |2011| ) (-1.5) or for a strip in the LMC using 
Herschel Science Demonstration Phase data by [Gordon et aL| 
( |2010 )). We thus decide to fix the effective emissivity index to 
1.5 and re-derive the temperature map. Fixing the emissivity 
index value prevents biases linked with 1 ) temperature-emissivity 
anti-correlation due to measurement uncertainties, 2) temperature 
mixing along the line of sight contributing to a bias in the derived 
"effective emissivity index" (see [Juvela & Ysard|20"T2l [Galametz] 
|et al.|2012[ for discussions on these issues). 

Error estimates - We perform Monte Carlo simulations in 
order to quantify the uncertainties on dust temperatures driven by 
flux errors of the resolved fluxes. Error maps are a combination 
in quadrature of the uncertainty maps derived during the data 
reduction (mapping, pointing) and uncertainties in the absolute 
calibration (15% for the SPIRE bands that are the most crucial for 
dust mass estimates), the latter being the dominant source of un- 
certainties. We generate 20 sets of modified constraints with fluxes 
randomly varying within their error bars and following a normal 
distribution around the nominal value and run the model for each 
set. Resolved uncertainties can thus derived using the standard 
deviations of each distribution. Since the absolute calibration is 
correlated for SPIRE bands, the three SPIRE measurements move 
consistently in each set. 



3.2 Realistic dust properties fitting 

We also apply the phenomenological SED fitting procedure of |Gal-| 
|liano et aL| (^2011) to our 3.6 to 500 /um data in order to derive the 
resolved total infrared (TIR) luminosities, dust properties and radi- 
ation field intensity estimates across the region. Their approach is 
similar to that of ^Draine et al.| ( [2007] ), except for the starlight inten- 
sity distribution. In both models, as suggested by |Dale et al.| ( |200l) , 
the distribution of starlight intensities per unit dust mass is approx- 
imated by a power-law (of index a). In the D raine et al.| ( [2007| ) 
models though, a significant fraction of the dust mass is assumed 
to be heated by starlight with a single intensity Vmin ("diffuse ISM" 
component). This model leads to better resolved fits (lower val- 
ues) for nearby spirals (see Aniano et al.|2012| ) but does not seem to 
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correctly reproduce LMC observation s ([Galliano et al.|2011| ). We 
choose the approach of |Galliano et al.| ( |201 1| ) in order to allow flex- 
ibility of the model at submm wavelengths and account for possible 
variations of the effective emissivity or the presence of colder dust. 

We consider the dust size distribution and composition to 
be uniform across the complex and that sources of IR emission 
are old stars and dust (PAHs, carbon grains and silicates). The 
stellar contribution to the MIR is modelled using a library of 
stellar spectra previously synthesised using the stellar evolution 
code PEGASE ( |Fioc & Rocca-Volmerange|p^97| ). The stellar 
population is assumed to have undergone an instantaneous burst 5 
Gyr ago (initial solar metallicity, Salpeter Initial Mass Function). 
We do not have a MIR spectrum to properly constrain the PAH 
composition. We thus choose to fix the ionised PAH-to-neutral 
PAH ratio (fpAH+), an additional parameter of the [Galliano et al.| 
( [201 l| l model, to 0.5 in order to limit the number of free parameters 
of the model. 

The free parameters of our model are thus: 

• the total mass of dust (M^ust), 

• the PAH-to-dust mass ratio (f^A//), 

• the index describing the fraction of dust exposed to a given 
intensity (of), 

• the minimum heating intensity (Umin)^ 

• the range of starlight intensities (AU), 

• the mass of old stars (scaling of our stellar spectra). 

We note that U=l corresponds to the intensity of the solar 
neighbourhood (2.2xl0~^ W m~^). For the rest of the paper, 
\Jmin+^^ will be referred as ''Umax''- The PAH fraction fpAH is 
normalised to the average Galactic value (4.6 %; [Li & Draine[ 
[2001) ). As pointed out in [Galliano et al.[ ( [20n] ), the mass of old 
stars is poorly constrained by the IRAC bands and is only used 
to subtract the stellar continuum to the MIR bands. No balance is 
made between this parameter and the starlight intensity. Errors on 
the physical parameters are derived using the same Monte Carlo 
technique described in Section 3.1. 

Amorphous Carbon to model carbon dust - Two diff'erent ap- 
proaches are compared in Galliano et al.[ ( [2011[ ) as far as the mod- 
elling of carbon dust is concerned: the "standard model" that uses 
the [Zubko et al.[ p004[ ) Galactic grain composition and graphite 
grains, or the 'AC model" that uses amorphous carbons (Zubko 
[et al.[[T99 6) in lieu of graphite. Carbon constitutes a major frac- 
tion of the grains in galaxies, especially in the circumstellar dust of 
carbon stars. However, the exact form of carbon dust is still rather 
unclear. Many SED models use graphite to describe the interstel- 
lar carbon dust as first suggested by[Mathis et al.[(|1977 ). Graphite, 
indeed, has well-known physical properties, which makes it easy 
to model, and is in good agreement with the observed (Galactic) 
extinction and polarization ( [Draine & Lee|1984| . However, obser- 
vational evidence challenges the graphite theory (at least the pure 
mono crystalline graphite theory), among which are 1 ) variations in 
the 2175 A profile not explained by changes in graphite grain size 
( [Draine & Malhotra,1993| ), 2) a graphite broadband emission fea- 
ture near 33 iim less strong (weaker and narrower) than expected 
in global models (e.g. Draine et al.'2007 ), 3} an erosion of carbon 
dust in shocks not reproducible in the case of graphite grains. In- 
deed, [WtyetaT] ([2^2]) or [Serra Diaz-Cano & Jones] ( |2008l ) show 
that graphite is not expected to survive as such in the ISM due to 
erosion and irradiation in shock waves and suggest that interstellar 
hydrogenated amorphous carbons could be the most probable form 



of carbon material, their erosion being more efficient than graphite 
in shocks. We will also see further in the paper that using the "stan- 
dard model" leads to dust masses inconsistent with the dust-to-gas 
mass ratio expected for the complex. The use of AC to model car- 
bon dust generally leads to a decrease of the mass of large grains 
compared to that required by the "standard model" to account for 
the same emission (see Section 4.5.1 for explanations). Consider- 
ing the average emissivity index of our region (y5~1.5), we favour 
the "AC" approach to model the resolved SEDs across the complex. 



4 RESULTS AND ANALYSIS 

In this section, we describe the various maps obtained directly or 
indirectly from our resolved SED models and analyse how their dis- 
tributions correlate. Figure [4] shows the maps of the parameters we 
derive (temperature and mean radiation field intensity, PAH frac- 
tion, star formation rate and dust mass surface density). The 8yL/m 
image is shown on the upper left panel for comparison. In Section 
4.4, we also investigate how submm colours evolve with the ISM 
conditions (heating sources, radiation field intensities). Finally, we 
analyse the dust surface density map and how it relates to the dis- 
tribution of the different gas phases in the complex, and more par- 
ticularly in the N159 region in Section 4.5. 



4.1 Temperatures and mean radiation field intensities 

We show the final cold temperature (i.e. temperature of the cold 
MBB) map obtained using our resolved two-MBB fitting technique 
in Fig. |4] (upper middle panel). The temperatures peak at the loca- 
tions of star-forming regions where dust is expected to be warmer. 
N160 shows the highest cold dust temperature, with a maximum of 
~40K while the median temperature of the N158-N159-N160 com- 
plex is 26.9-^^ K (28.2 if we restrict the analysis to ISM elements 
with a 3-cr detection in the SPIRE bands). Based on observations 
made with the TopHat telescope combined with DIRBE data (100, 
140, and 240 //m, 42'), [Aguirre et al.[ ( [2003 ) derived temperatures of 
25=^i-^K and 162^^^}^ for the LMC and 30Dor, using 1.33^^ ^^ 
and 1.5-^^^, respectively. From Spitzer data (thus no data above 160 
yL/m) and the same emissivity values, ,, Bernard et a"r] ([2008 ) obtained 
a colder temperature of 21.4 K for the LMC and 23 K for the 30 
Dor region. Our values are thus close to the values derived for the 
30Dor region by [Aguirre et al.[ ( [?003[ ). We note however that diff'er- 
ent methods were used to derive the temperatures, making a direct 
comparison difficult. The N159S region is finally the coldest region 
of the whole complex, with a temperature ~22K. This supports the 
argument that massive stars may still be at infancy and deeply em- 
bedded in the circumstellar material in this particular region. 

In order to investigate more precisely the distribution of ra- 
diation field intensities across the region, we derive a map of the 
mass-weighted mean starlight heating intensity <U> ( [Draine & Li[ 
[2007[ ) that develops into: 

AU 
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with AU the range of starlight intensities, U^m the minimum 
heating intensity and a the index describing the fraction of dust 
exposed to a given intensity, all three obtained from the resolved 
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Figure 4. Results of our resolved SED fittings: Upper left: 8 //m map convolved to the SPIRE 500 //m resolution (units of MJy sr~^). Upper middle: Cold 
dust temperature (units in Kelvins). Upper right: Mean starlight heating intensity (U=l corresponds to the intensity of the solar neighbourhood). Lower left: 
PAH-to-total dust mass ratio normalised to the Galactic value (4.6%). Lower middle: Star formation rate (units in M© kpc~^ yr~^ - log scale). The white cross 
indicates the N159S region. Lower right: Dust surface density (in M© pc~^). For all images, North is up, east is left and coordinates (RA, Dec) are given in 
equinox J2000. 



|Galliano et al.| p011| ) modelling procedure. U= 1 corresponds to 
the intensity of the solar neighbourhood (2.2xl0~^ W m~^). 

We show the mean starlight heating intensity map in Fig. |4] 
(upper right panel). Since dust grains are heated by the interstel- 
lar radiation field, we obtain a similar distribution between dust 
temperatures and radiation field intensities. Maxima in <U> ap- 
pear close to Hii regions where the starlight heating the dust is 
expected to be much more intense than the overall starlight illumi- 
nating the complex. The averaged intensity in our ISM elements is 
<U>=8.6. This average increases to 11.7 when we restrict the cal- 
culation for resolved ISM elements that have a 3-cr detection in the 



SPIRE bands. We note thatl Meixner et al.|p010| ) derived an average 
value <U>=7.6 for a strip across the LMC (so including more qui- 
escent regions) using the same model as that applied in our study. 
Finally, N159S shows the lowest <U> (~1.8) of the complex. This 
is consistent with the cold temperatures we observe in this region. 
This region thus appears to be shielded from the strong starlight 
heating the N159 Hii region. 



4.2 PAH fraction 

Figure |4] (lower left panel) shows how the PAH-to-total dust mass 
ratio fpAH is distributed across the star-forming complex. We re- 
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mind the reader that in the modelUng, fpAH is normaUsed to the 
Galactic value fpA//,Mw^=4.6%. The average value of fpAu across 
the region is driven by the ISM outside the H ii regions themselves 
(median fpA//=0.7-^^). However, we observe variations in the PAH 
fraction. Peaks of the fpAH (in red) are located in low-surface bright- 
nesses regions and are probably artefacts due to the very low dust 
mass of those regions. We observe a strong decrease of fpAu in 
the three star-forming regions as well as low PAH fractions in two 
compact knots located at the south east side of N160. All these re- 
gions correspond to peaks in the radiation field intensities (Fig. |4] 
<U> panel) as well as strong radio emission (as we will see later 
in Section 5). This depression of the PAH fraction by intense ra- 
diation fields was also studied in detail by |Sandstrom et al.1 ( |2010| ) 
in the SMC. These results are consistent with the fact that PAHs 
are known to be destroyed in strong Hii regions ( [Madden et al.| 
|2006[ [Lebouteiller et al.|200 7V PAHs are also smaller on average 
in low-metallicity galaxies, which could also favour their destruc- 
tion ( [Sandstrom et al.'2()12a). Low-metallicity environments often 
show harder interstellar radiation fields extending over larger size 
scales. This eff'ect is attributed to the hardness of the radiation field 
and/or to a lower dust attenuation in low-metallicity ISM which 
consequently leads to a larger mean free path length of the ionis- 
ing photons. Moreover, the fragmented structure and dumpiness 
of the ISM in low-metallicity galaxies allows UV light to penetrate 
deeper in metal-poor molecular clouds, favouring grain destruction. 
We note that should the radiation field be harder than the interstellar 
radiation field of Mathis et al. ( 1983 ) used in our model, the PAH 
fraction we estimate could be over-estimated. 

|Dwek|p005l ) and [Galliano et al.| ( [2008 l) have suggested that 
the low fraction of PAHs in metal-poor environments could also 
be due to a delayed injection of carbon dust into the ISM by AGB 
stars in the final phase of their evolution. Using near-IR JHK pho- 
tometry in the N159 nebula, [Meynadier et al.[ ( [2004] ) distinguished 
a young (~3 Myr) and an old (~1-10 Gyr) stellar population, even 
if they could not exclude that both populations might be spatially 
unrelated. Based on the absence of red supergiants in N159 (stars 
that should normally be observed in the cluster if the present stellar 
population was at least 10 Myr old), [Jones et al.| ( [2005[ ) concluded 
that most of the star formation must have taken place recently and 
that N159 might not be older than 1-2 Myr. A sequential cluster 
formation from N160 to N159S has been observed by [Nakajima[ 
[et al.[ p005[ ), using observations of Herbig Ae/Be clusters. This 
supports a slightly longer star formation activity for the N160 
region (~ 10-30 Myr). Both regions are nevertheless quite young 
compared to the time necessary for carbon dust to be injected in the 
ISM by the stars and thus participate in the formation of PAHs in 
the clusters. We finally note that recent studies suggest than the rate 
of AGB -produced PAHs compared to the rate of PAHs formed in 
the ISM may also change with metallicity ( [Sandstrom et al.[2012a[ ). 

In conclusion, although the obvious anti-correlation between 
fpAH and the radiation field favours the destruction of PAHs as a 
straight forward explanation for the low PAH fraction we observe, 
the young age of the clusters could be an additional explanation for 
the absence of PAH in our three H ii regions. 

4.3 Star Formation rate 

Since the total bolometric IR luminosity traces the emission of stel- 
lar populations enshrouded in dust cocoons, this quantity is thus 
very often used to quantify the star formation obscured by dust 
( [Perez-Gonzalez et al.[[20ngl [Kennicutt et"aLl[2009l ). We first ob- 



tain the total IR luminosity map of the complex by integrating the 
resolved SEDs in a v-fy space: 

XI 100 //m 
Ly dv (3) 

In the SED fitting technique we use, we model the stellar con- 
tribution to short wavelengths as a separate component. This contri- 
bution is subtracted before performing the integration to only take 
the thermal emission of dust into account in the calculation. We 
use the IR-based calibration of [KennicHttl ( [1998) to derive a star 
formation rate map of the region: 



Mo jr-i 5.8 X lO^Lo 

[Jones et aL| ( [2005| ) estimated the integrated luminosity of the 
N159 complex (integration of the SED up to lOOyum combining 
Spitzer and KAO data). They found that the integrated luminosity 
was consistent with the observed radio emission, assuming that 
the Hii region was a cluster with a normal Initial Mass Function 
(IMF). As mentioned in [Indebetouw et al.[ ( [2008[ ), Eq. |4] was 
calibrated for large extragalactic sources. It may be valid for the 
full star-forming regions such as N159, but the SFR calibration 
may not be accurate when applied to some of our resolved ISM 
elements, namely those not properly covering the stellar IMF. 

We show the SFR map of the region in units of M© yr~^ kpc~^ 
(log scale) in Fig. [4] (lower middle panel). We find a median SFR 
of 0.06 Mq yr~^ kpc~^ across our modeled region. This average 
increases to 0.17 M© yr~^ kpc~^ when we restrict the calcula- 
tion for resolved ISM elements that have a 3-cr detection in the 
SPIRE bands. This is consistent with the average SFR obtained 
by In debetouw et al.[ ( [2008j using the H2 surface densities and 
the Schmitt-Kennicutt law (0.14 Mq yr'^ kpc'^ for N159/N160 
and 0.11 Mq yr~^ kpc~^ for the southern molecular ridge to whom 
the complex belongs to). They find that their SFR estimate for 
N159/N160 corresponds to that derived using the Calzetti ( 2007]) 
formula (i.e. a calibration from Ha and 24 jim luminosities) within 
uncertainties. We note that the SFR we derive varies significantly 
across the region. Peaks of SFR correspond to our three star- 
forming regions, with the N160 centre showing the highest star 
formation rate of the complex. We indicate the N159S region by 
a white cross in Fig.|4](SFR panel). The region does not show mas- 
sive on-going star formation. 

4.4 Submillimeter colours 

[Bendo et al.| ( |20T2l l studied Herschel colours (i.e. PACS or SPIRE 
surface brightness ratios) in three nearby galaxies (M81, M83 and 
NGC 2403) and found that the 70/160 colour was correlated with 
tracers of the star formation activity (thus that warm dust is pri- 
marily heated by the young stellar populations) while SPIRE ra- 
tios (250/350 or 350/500) were more strongly correlated with the 
galacto-centric radius of the galaxies. They concluded that cold 
dust could be primarily heated by the total stellar population (young 
stars + old stars). [Boquien et al.[ ( [20Tm and'G alametz et al.[ ( [20T0] l 
studied similar correlations in M33 and NGC 6822 respectively and 
showed that these ratios were also strongly correlated with the star 
formation, traced by the 24 yt/m emission for instance. In Section 
4.4.1, we thus study the dependence of FIR colours on various pa- 
rameters, including the mean radiation field intensities <U> and 
the PAH fraction just derived from our modelling. Moreover, if 
FIR colours are often used as good proxies for dust temperatures. 
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Figure 5. Comparison of PACS and/or SPIRE surface brightness ratios to the 3.6 fim surface brightnesses, the 24 jum surface brightnesses, the mean stellar 
density <U> or the PAH fraction fpAH- Black lines show the best fit lines and we provide the slope and the RMSE at the bottom of each panel. 



colour-colour relations can, on the other hand, provide additional 
information about the dust emissivity wavelength dependence (|Na-| 



Igata et al.|2002||Hibi et al.|2006||Hrrashita & Ichikawa|2009| ). In 
Section 4.4.2, we analyse SPIRE colours as possible indicators of 
the cold dust emitting properties. 

4.4.1 Dependencies 

In order to investigate the drivers of submm colours variations in 
the N158-N159-N160 complex, we compare in Fig. |5] the PACS 
and SPIRE surface brightness ratios 100/160, 160/250, 250/350 and 
350/500 with : 

1) the 3.6 fim surface brightness: dominated by emission from 
the old stellar populations, it is known as a good proxy of the stellar 
mass (see |Zhu et al.|2010l|Eskew et al.i2012 for recent studies), 

2) the 24 iim hot dust emission, often used as a good calibrator 
for star formation ( |Calzetti|2007[|Rieke et al.|2009| among others), 

3) the mean starlight heating intensity <U>. Values of U have 
been normalised to that of the Milky Way), 

4) the PAH fraction fpA//. 

For this section as well as the following, we restrict ourselves 
to pixels with a 3-cr detection in the three SPIRE bands. We note 
that part of the N159S region matches this signal-to-noise criteria. 
In Fig. [5] we compute the linear regressions from each comparison 
using an ordinary Least-Squares regression (function linfit.pro of 



the Astrolibrary of IDL). On each plot, we provide the slope of the 
relation as well as the root-mean- square error defined as: 



RMSE : 



^{Oi-Eif 



(5) 



with the observed ratios (or temperature colours), E/ the ratios 
estimated from our fitting relations and n the number of ISM 
elements. 

We observe correlations in all cases, with a very tight rela- 
tion between the lOO-to-160 and 160-to-250 flux density ratios and 
the mean starlight intensity <U> (lower RMSE for each line). We 
know that dust grains are heated by the interstellar radiation field 
and re-emit this energy in the infrared regime so correlation be- 
tween FIR colours and the radiation field intensity are expected. 
FIR colours are often used as good proxies of the average dust 
temperatures, even if FIR ratios using shorter wavelengths than 
100 jim (70-to-lOO ratio for instance) can be potentially biased by 
emission from stochastically heated dust grains ( [Dale et al.|2001[ 
[Bernard et al.pOOS] ). Contrary to ratios at shorter wavelengths that 
show stronger correlations with the 24 jim surface brightnesses, the 
SPIRE-only ratios follow the 3.6 and 24 fim surface brightnesses 
equally well. We finally note that the N159S region shows the low- 
est submm surface brightness ratios of the whole complex, with 
surface brightness ratios of -0.75 for 100/160, -1.58 for 160/250, 
-2 for 250/350 and -2.35 for 350/500. 
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Figure 6. Colour-colour plot of the star-forming complex: 350/500 colours are plotted against the 250/350 colour. Squares represent the resolved ISM 
elements of the star-forming complex (pixel size = lA" i.e. ~4pc for the LMC) coded with temperature (temperatures taken from the map presented in Fig. [4}. 
We separate the resolved elements in 4 panels (T<26, 26<T<28, 28<T<30 and T>30). For each panel, elements are binned 3-by-3 for clarity. The lines show 
the colours obtained from a single temperature modified blackbody with an emissivity index (SsingieBB of 2, 1.5 and 1, from cold temperatures to the left to hot 
temperature to the right. 



4.4.2 SPIRE colour-colour diagrams 

Figure [6] shows a SPIRE colour-colour diagram of the star-forming 
complex. We plot the 350-to-500 fim flux density ratios (hereafter 
350/500 colour) of our resolved elements as a function of the 250- 
to-350 fim flux density ratios (hereafter 250/350 colour). In order 
to investigate variations in the submm colours with temperature, 
we separate the points in 4 panels (T<26, 26<T<28, 28<T<30 and 
T>30) and code them with temperature, from cold in dark purple 
to hot in red colours. We use the resolved temperatures derived 
using our two-modified blackbody model. The lines show mod- 
elled colours obtained from a single temperature modified black- 
body with an emissivity index psingieBB of 2, 1.5 or 1. 

We observe that SPIRE colours are well correlated with each 
other. Resolved elements with T>30 (fourth panel) have, on aver- 
age, slightly lower 350/500 colours for a given 250/350 colour, so 
a "flatter" submm slope. This trend agrees with the anti-correlation 
observed between the emissivity index and the temperature (Du- 
|pac et aL|2003l|Shetty et aH 2009 ) but could also result from tem- 
perature mixing eflTects along the line of sight in these elements. 
Given our error bars, we, however, consider this trend as marginal. 
We note that the ISM elements showing very low 250/350 colours 
(black squares) belong to N159S. 

As suggested by Fig. [6] a single modified blackbody model 
with PsingieBB-'^ cannot reproduce the colours we observe. In fact, 
many resolved elements reside around the p-\ model, thus at the 
very low end of the p values quoted in the literature {\<p<2.5). 



This cautions the use of single-temperature fits to derive the intrin- 
sic emissivity of the dust grains, those models providing in reality 
an "eflfective" emissivity index resulting from the combination of 
various dust populations at diflferent temperatures. 

4.5 Dust and gas masses across the complex 

4.5.1 Dust masses 

We show the dust surface density map l^dust obtained with the "AC 
model" in Fig. |4] (lower right panel). The distribution of the dust 
surface density follows the distribution of the FIR/submm emis- 
sion. We particularly point at the elongated structure of N 159 cor- 
related with the SPIRE and 870 yt/m emission. The total dust mass 
of the region we modeled is -2.1x10^ M©, with 8.5x10^ M© within 
resolved elements fulfilling the 3-cr SPIRE-band detectior]^ A sig- 
nificant amount of dust is present in N159S, a region whose emis- 
sion starts to be detected at 160 jim with PACS (Fig. [2|, thus a 

^ We note that estimating this total dust mass on a resolved scale (resolved 
masses added together) allows us to probe the various dust components 
in ISM elements with low-dust mass surface densities as well as in dense 
regions. Contrary to dust mass estimates derived from a global SED mod- 
elling of the region, we can avoid potential biases (namely a dust mass un- 
derestimation) caused by a poorer resolution (c.f. discussion about biases 
caused by non-linearities in SED models in |Galliano et al.|2011| . 
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significant reservoir of cold dust. We refer to Section 5 for an indi- 
vidual SED modelling of N158, N159, N159S and N160 that allows 
us to derive individual dust masses for these Hii regions. 

The resolved dust masses we estimate are 2.8 times higher if 
we use the "standard model" of |Galliano et al.|p011| ), i.e. graphite 
grains to model carbon dust. This is consistent with the statistical 
distribution of the resolved dust mass ratio between the "AC model" 
and the "standard model" they find in a strip across the LMC (me- 
dian of the distribution at ~2.6). In fact, the opacity profile of AC 
is higher than graphite for A>5fim. This implies that for a fixed 
starlight intensity, AC will have a lower temperature than graphite, 
thus lead to more mass. However, AC has a flatter submm slope 
(more emissivity in the submm regim^, so we are able to fit the 
same observational constraints with slightly hotter grains but less 
mass. Indeed, the "AC model" invokes higher starlight intensities 
(median of <U> across the complex: 8.6 in the "AC case", 4.1 in 
the "standard model"), and thus less dust than graphite to account 
for the same submm emission. We finally note that the errors on 
resolved dust masses are of the order of ~50%. This will introduce 
non-negligible uncertainties in our gas-to-dust mass ratio (G/D) es- 
timates (c.f. Section 4.5.4). 



4.5.2 Comparison with the Hi distribution 

High-resolution maps of the atomic gas reservoir have been built 
for the complete LMC. The H i data cube used in this section (from 
|Kim et al.|2003| ) combines data from ATCA, the Australian Tele- 
scope Compac t Array ( Kim et al.||2003|) and the Parkes single- 
dish telescope ( [Staveley-Smith et al .'2003'). The individual channel 
maps have a velocity resolution of 1.649 km s~^ and the Hi map 
was obtained integrating the cube over the range 190 < Vhei < 386 
with a resolution of T. We refer to Kim et al. (2003 1 



9x10' 



km s 



and 



[Bernard et al. ( 2008 ) for further details on the LMC integrated map. 

Figure |7] (top) shows the Hi map of the N158-N159-N160 
complex, with 870 jim contours overlaid. The distribution of H i 
does not match that of the emission in the FIR (both Spitzer and 
Herschel) bands across the region. A clear off'set is for instance ob- 
served between N159 and the closest H i peak located further north 
(~1.5'). A closer look at the regions indicates that the peaks of the 
870 jim emission of the three regions N158, N160 and N159, as 
well as the PACS and SPIRE peaks, reside in fact in H i holes. This 
could be a signature of the H i to H2 transition already observed in 
LMC or Galactic clouds ( |Wong et al.|2009 [ |Lee et al.|20"T2l Tat- 
ton et al. in prep, among others), with Hi being converted into H2 
either by thermal or gravitational instabilities or by shock compres- 
sions. Contrary to the other regions of the complex, the southern H i 
peak does spatially coincide with a peak in the 8 jim emission (see 
Fig. [3]bottom panel for comparison between the 8 and the 870 fim 
emission in N159S). This peak is located ~ 30 pc east of N159S. 



4.5.3 CO in N159 and N159S 

In this section, we focus more particularly in the N159 region to 
study how the dust budget and the dust properties are related to the 
molecular gas. Several CO observations of the N 159 nebula have 
been carried out so far. Johansson et al. (1998 ) observed ^^C0(J=1- 
0) with the SEST (Swedish-ESO Submillimeter Telescope). They 
obtained maps with a resolution of 45'', corresponding to ~1 Ipc for 



^ c.f. further discussions about opacities in 
A 



Galliano et al, 



j201l) , Appendix 
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Figure 7. Top panel: The N158-N159-N160 complex observe d in H i (in 
units of atom cm~^). The Hi map is a combination from ATCA jKim et al"] 
[2003) and Parkes jStaveley- Smith et al.|2003) observations. 870 jim con- 
tours are overlaid (0.02, 0.1, 0.25 and 0.4 Jy beam~^). Bottom panel: ASTE 



CO(3-2) observation of N 159 (in units of K km s~ 
top. For both images, North is up, east is left. 



). Contours: same than 



Nl 59. |Bolatto"erari ( [2OOOI 1 also mapped the N158 / N159 / N160 
complex in four transitions: ^^CO{\-0\ i2cO(2-l), ^^C0{A-3\ 
[Ci](3Pi ^ ^Vo) also with SEST and the 1.7m AST/RO (Antartic 
Submillimeter Telescope and Remote Observatory). A spatially re- 
solved survey of the LMC was also performed with the NANTEN 
4-m telescope in the 2.6 mm carbon monoxide emission with a 40 
pc resolution by |Fukui et al.| ( |1999l|2008| ), permitting a study of the 
giant molecular clouds at global scales. Several observations of the 
CO(4-3) transition and of the CO(3-2) transition have been carried 
out in selected regions of the LMC. These transitions have been ob- 
served in the N158 / N159 / N160 co mplex by|Mizuno et aLlpOTO] ) 
(i2cO(J=4-3) and i3cO(J=3-2) and Minamidani et al.|(|2008|) wi"^ 



NANTEN2 (beam size at 490 GHz: 18'') and ASTE (Atacama Sub- 
millimeter Telescope Experiment, beam size at 350 GHz: 22''-23''). 
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Figure 8. Left: Map of the gas-to-dust mass ratio of the N159 region derived using Xco=^ x 10^^ H2 cm~^ (K km s~^)~^ to convert CO observations in H2 
mass. 870 //m contours are overlaid (0.02, 0.1, 0.25 and 0.4 Jy beam~^). Right: Scatter in the logarithm of G/D across the region as a function of the X cp facto r 
used to convert CO observations in H2 mass (in units of M© pc~^(K km s~^)~^). We use the minimization technique described by jSandstrom et al.|j2012b) . 
The scatter is minimized when Xco = 8.6 Mq pc~^ (K km s~^)~^ . 



Finally, the Magellanic Mopra Assessment (MAGMA) survey has 
recently completed a high angular resolution ^^CO (1-0) mapping 
survey of giant molecular clouds across the LMC. The beam size 
of the instrument is 33'', which corresponds to a spatial resolution 
of ~8 pc at the distance of the LMC. 

The two giant molecular clouds of N159 (N159E and N159W) 
as well as its southern region (N159S) are thus already mapped at 
high resolution in many transitions. Figure^ shows the N159 re- 
gion observed by ASTE CO(3-2) with 870 //m contours overlaid. 
N159W is the strongest region and the eastern region breaks into 3 
independent peaks. Compared to the other giant molecular clouds, 
the N159S does not show any heating source and little star forma- 
tion activity. [Bolatto et aL| ( |2000| ) also note the detection of diffuse 
low level emission in C0(2-l). 

N159E and N159W are associated with embedded star clus- 
ters with temperatures that were determined to be 70-80K while 
N159 S shows a nearly uniform molecular temperature estimated 
to be ~30K by |Mizuno et al.| ( [201 ). They also note that the molec- 
ular peak of the N159E region is showing an elongated structure 
that is very well resolved at 870 jim with LAB OCA. A slight shift 
appears between the peak of ^2CO(4-3) in N159W and the 8 and 24 
yL/m peaks. The bright peak of the LABOCA emission in N159W 
does not seem to show any shift with the ASTE ^^CO(3-2) map. 
They also compare the spatial distribution of CO with Spitzer bands 
and found a fairly good correlation between the maps. As men- 
tioned before, a peak in the 8 yum emission, situated on the east 
side of N159S, corresponds to a peak in the Hi distribution while 
a peak in the 870 jum emission (N159S) corresponds to a peak in 
the CO map. [Fukui et aL] p009| ) already noted the offset between 
the Hi and CO peak for N159. They showed that the Hi and CO 
distributions correlate well on a 40-100 pc scale and argued that the 
Hi envelopes are gravitationally bound by giant molecular clouds. 
They also show that the correlation between the Hi and CO dis- 
tribution breaks down on a more resolved scale. They interpret it 
as an illustration of the conversion of warmer, low-opacity Hi to 



colder high-opacity H i from which H2 could form, as suggested in 
|Ott et aL| ( [2058l ) and |Wong et aLl ( [2009l ). 

4.5.4 Gas-to-dust mass ratios in N159 

Assuming a close mixing of dust and gas and that the gas-to-dust 
mass ratio does not depend on whether we study a predominantly 
atomic or molecular phase, we can relate the dust mass surface den- 
sity in each resolved element of N159 to the gas reservoir such as: 

GID = ^"'^^colco 

^dust 

with Xco the conversion factor in units of M© pc~^ (K km 
s~^)~^ and E/// and ILdust the Hi and dust surface densities in 

Mq pC-2. 

We first need to convert the ASTE CO(3-2) observation of 
the region N159 (shown in Fig. [7] bottom panel; private commu- 
nication) to a CO(l-O) map. A mean CO (J = 3-2)/CO (J = 1-0) 
ratio (R3_2/i-o) of 0.6 is usually derived in galaxy disks ( |Warren| 
et al.|2010| but R 3_2/i-o can probe larger ranges (0.2<R3_2/i_o<2; 
Meier et al.|200T||Tosaki et aL|2007[rMao et al.|2010| [Warren et"aLl 
2010| ). [Minamidani et aL] (^2008) investigate this ratio in 10 CO 
clumps of N 159. This Hii region being a hot and low-metallicity 
environment, they find values ranging between 0.7 and 1.4 (so 
higher than the median 0.6 value), with R3_2/i-o=0.7 in the core 
of N159S and 1.1 in the cores of N159W and N159E. No clear 
North-South trend is observed however. We use the average 
R3_2/i-o value of the 10 clumps (~0.9) to convert the CO(3-2) into 
CO(l-O). From the H i map and the derived CO(l-O) map, we can: 

Derive a gas-to-dust mass ratio map - We use a constant Xco 
factor to convert the CO intensity in molecular hydrogen mass 
surface density. |Fukui et al.| (^009 ) has performed a survey of 
molecular clouds in the LMC and estimate Xco to be ~7 x 10^^ H2 
cm~^ (K km s~^)~^ , so three times the Galactic value of ~2.3 x 10^^ 
H2 cm~^ (K km s~^)~^ . We use the dust mass map obtained with the 
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"AC model" as our dust reference map and convolve the CO and 
dust surface density maps to the resolution of the Hi map using 
Gaussian kernels. Figure [8] (left) shows the gas(Hi+CO)-to-dust 
mass ratio map obtained with those hypothesis, with a mean 
gas-to-dust mass ratio (G/D) value of 384 (with a non-negligible 
uncertainty given the 50% errors on dust masses). As a reference, 
the Galactic value of G/D is -157 ( |Zubko et aLl[2004l ) so the 
value we derived is 2.4 times the solar value. The total Hi mass 
contained in the colored ISM elements of Fig. [8]is 4.3 X 10^ Mo 
while the molecular gas mass is 1.1 x 10^ M© (thus 2.6 times 
higher than the atomic mass). The total dust mass in these ISM 
elements is 4.3 x 10^ M©, leading to a globally derived average of 
356, consistent with the value derived on a resolved basis in spite 
of the local variations. We note that given the higher contribution 
of the H2 to the total gas budget, using a slightly different (but 
constant) Xco factor would simply translate into a crude scaling 
of Fig. [8] to higher or lower gas-to-dust mass ratio ranges but 
would barely affect the qualitative distribution of Fig. [8] (left). 
Finally, dust masses are 2.8 times higher when graphite is used 
to model carbon dust ("graphite model"), leading to G/D smaller 

[Franco & Cox|(p^86l), 



by the same factor (G/D~127). Studies of 



[Lisenfeld & Ferrara| { [1998| ) or [James et al 



([2002)), among others, 
show that the G/D is proportional to the metallicity of the galaxy. 
The G/D obtained when dust masses are obtained with graphite in 
the SED modelling are thus too low compared to those expected in 
metal-poor environments like the LMC. 

Try to constrain the Xco factor - If we assume now that G/D 
is constant across the region, we can use the technique proposed 
by Sandstrom et al. ( 2012b ) in order to derive the Xco factor that 
minimises the scatter of the logarithm of the G/D across the region. 
This technique minimises the effects of outliers on the measured 
scatter. As suggested in San dstrom et al.[ ( [2012b| , we tabulate the 
Xco values from 0.1 to 100 M© pc'^ (K km s'^)"^ (thus from 6.2 x 
10^^ to 6.2 X 10^^ H2 cm"2 (K km s"0"0 and estimate the scatter 
of the logarithm of the G/D for each value of Xco. The scatter for 
a given Xco value is estimated using the IDL biw eight -mean.pro 
function. We finally take the Xco value that minimises the scatter 
(Figjs]- right). The minimization occurs for Xco =8. 6 pc~^ (K 
km s-^-i (= 5.4 X 10^0 H2 cm-^ (K km s-^)-^, thus 2.3 times the 
Galactic conversion factor and 1.3 times lower than that derived by 
jFukui et aL ( 2009|). Th is is slightly higher than the factor derived 
by 'Pine da et al.[ ( [2009| ) using CO observations only (and the virial 
mass/density profile of CO clumps) for the m olecular ridge or the 
dust-based value derived by Leroy et al. ( |201l| ) (~4 x 10^^ H2 cm~^ 
(K km s"0"0 for the LMC using Spitzer/SAGE data only. 

We caution the fact that [Sandstrom et al.[ ( [2012b| ) apply the 
minimization method to unresolved molecular clouds. Xco is, in 
that case, a statistical quantity involving a cloud ensemble, namely 
a measure of the covering factor of CO clouds in more extended 
molecular structures. This method favors CO-bright gas and ig- 
nores CO-dark gas, which could naturally lead to Xco values close 
to those of the Solar Neighbourhood. This technique may not be 
applicable to the individual molecular clouds or clouds complexes 
we resolve in the LMC. Further studies on the Xco factor and po- 
tential CO-dark gas reservoirs using dust masses as a tracer of the 
gas masses will be performed in a future study. 



5 THERMAL EMISSION AT 870 fim 

As mentioned previously, the LAB OCA data reduction might re- 
move faint diffuse emission in the outskirts of the complex. In this 



section, we select 24 regions across the complex where the 870 jum 
emission is bright enough not to be affected by the filtering steps 
and analyse the thermal dust emission at 870 yum in those regions. 

5.1 Selected regions and correction of radio contamination 

We select individual regions across the complex that span a wide 
variety of environments, from intensely star-forming regions to a 
more quiescent ISM to study the detailed dust properties with the 
empirical SED model of [Galliano erari ( [20TT) . Figure [9] (top left) 
shows those regions. The circles indicate the photometric apertures 
we use. We detail the characteristics (centers and sizes) of those 
apertures in Table [T] along with the Herschel and LAB OCA 870 
jjm flux densities of the individual regions. Herschel errors are es- 
timated using a combination in quadrature of calibration uncertain- 
ties and measurement uncertainties (background effects, choice of 
the aperture). LAB OCA errors are quantified from the rms maps 
produced by the reduction pipeline. 

The 870 fim flux densities we quote in Table [T] contain 
non-dust contribution, especially free-free radio emission. 4.8 
(resolution 30'') and 8.6 GHz (resolution 15'') images of the entire 
LMC have been obtained by Dicke l et al. (2005 ) using ATCA, 
which allow us to quantify the radio contamination. These obser- 
vations reveal a strong radio emission in 30 Dor (peak of the radio 
emission of the entire LMC) as well as in the 3 bright star-forming 
regions N158, N159 and N160 (see Fig. [9] top right). They also 
indicate radio emission in two compact knots in the south east 
of N160, spatially correlated with peaks in the IR emission (as 
well as peaks in the SFR or <U> maps) and corresponding to 
local decreases of the PAH fraction (see Fig. |4|. We convolve the 
ATCA maps to the SPIRE 500 fim resolution using a Gaussian 
kernel and estimate the 4.8 and 8.6 GHz fluxes (at 6.25 and 3.5 cm 
respectively) within our apertures for the 7 regions showing radio 
emission, namely region 1,1a, lb, Ic, 2, 3 and 5. With these two 
constraints, we perform a regression analysis (Lradio ^ v~^ and 
estimate a free-free contamination of ~9.6 %, 9.0 % and 13 % of 
the 870 jum fluxes in the regions 2, 3, 5 respectively and 10.4 % in 
region 1 (N159), with individual contaminations of 8.0 %, 13 % 
and 5.6 % in the subregions la, lb and Ic respectively. 



5.2 Results 

We first apply our two-MBB model (J3c=l.5) to derive the cold dust 
temperature of our selected regions using the integrated flux den- 
sities tabulated in Table [T]but correcting the 870 jum flux for radio 
contamination. We also model the selected regions using the "AC" 
model of [Galliano et al.| ( [2011 1) described in Section 3. For the 7 
regions where radio emission is detected, we kept the 870 fim flux 
densities of Table[T|(i.e. flux densities including non-dust contami- 
nation) but include the radio component in the modelling procedure 
to allow the visualization of this contribution in the SEDs of those 
regions. This component will be constrained by our 2 ATCA data 
points. As mentioned in the previous section, the free-free contam- 
ination is inferior to 10% in the Hii region and negligible outside 
their cores. 

Figure[9] (bottom panel) shows the SEDs of regions detected in 
radio (and modelled with a free-free component) while Figure [T0| 
shows the SED obtained in the other regions. We observe varia- 
tions in the SED shapes from region to region. The SEDs of the 
brightest H 11 regions (regions 1, 2, 3, 5) show 'flat' profiles indicat- 
ing a warmer dust temperature range. The main differences in the 
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Figure 9. Top left panel: Selected regions across the complex overlaid on the 870 yum image. Circles also indicate the size of the photometric apertures. Top 
right panel: ATCA radio map at 4.8 GHz (see Dickel et al.|2005| for the complete map of the LMC). Bottom panel: Individual SEDs of the bright sources 
using the Galliano et al. ( 201 1) "AC" SED modelling technique in addition with a radio component. The SEDs models shown in this figure were thus produced 
using amorphous carbon to model carbon dust. Filled circles show Spitzer, Herschel, LABOCA and ATCA radio fluxes. Errors are smaller than the symbols. 
The integrated SEDs appear with thick black lines. The red lines underline the thermal dust contribution. The dash-dotted lines indicate the free-free emission. 
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Figure 10. Individual SEDs in different regions of the diffuse ISM using the |Galliano et al.|j2011) "AC" SED modelling technique. The SEDs models shown 
in this figure were thus produced using amorphous carbon to model carbon dust. Filled circles show Spitzer, Herschel and LAB OCA fluxes. Errors are smaller 
than the symbols. 



SED shape occur in the mid-infrared regime where the 24-to-70 
slope strongly varies from one region to another, with higher 24/70 
colour in H ii regions and lower ratios in quiescent ISM (regions 12, 
19 or 21 for instance). Region 6 (N159S) and region 21 associated 
with the same molecular cloud have SEDs strongly peaked due to 
the absence of star formation in the region. 

Table |2] summarises the most important parameters obtained 
from the two fitting procedures: the cold temperature T^iust^ the dust 
mass Mdust, the PAH fraction fpAu and the three parameters charac- 
terising the radiation fields (the minimum and maximum intensities 



of the radiation field \]min and \5max and a, the index describing the 
fraction of dust exposed to a given intensity). We also add the dust 
masses obtained when a "graphite model" is applied to the data, for 
comparison. As mentioned previously, diffuse ISM shows dust tem- 
peratures of ~26K on average while N158, N159 and N160 show 
warmer dust temperatures (>30K). There is thus a clear evolution 
of the dust temperature range between bright star-forming regions 
like N159 (region 1), intermediate star-forming regions like region 
3, 17, 16 to very quiescent regions like region 10 or 21. The dust 
masses are systematically lower when using amorphous carbon in 
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lieu of graphite to model carbon dust, 2.8 times lower on average. 
Indeed, the "amorphous carbon" model requires less cold dust to 
account for the same submm emission. This will have significant 
consequences on the study of the gas-to-dust mass ratio in the star- 
forming complex (see Section 6 for discussion). We finally note that 
the values derived using graphite are of the same order than those 
derived by |Rantakyro et al.| ( [2005 l) usins SIMBA data at 1 .2mm. 

The a values for the diff'use regions very often reach the limit 
of 2.5 fixed in the model for this parameter. This is consistent with 
the assumption that in the diff'use medium the heating intensity 
falls as r~^ away from the heating source so dU/dr evolves as r~^ . 
In a uniform medium where dMj^^^/ dr evolves as r^, dMj^^^/dU is 
thus proportional to U~^^, so a-2.5 ( [Dale et al.|200l|). To perform 
a global SED modelling of nearby spiral galaxies, [Praine et alT] 
P007| ) and [Dale et al.| ( [20T2] ) advise fixing the maximum heating 
intensity to 10^, to limit the number of free parameters (fixed to 10^ 
injAniano et al.|2012|). Given the number of data points we have to 
perform the modelling, we decide not to fix the maximum heating 
intensity of the radiation field in our study to investigate its 
possible variations from one region to another. Our Vmax values are 
not allowed to be higher than 10^ during the modelling. Only two 
regions reach this maximum limit. N158, N159 and N160 show 
maximum intensities of ~1 0^-10^ while diff'use ISM shows lower 
maximum intensities. As suggested by their peaked SED shape, the 
two south regions N159S and Region 23 have a very narrow (and 
low) range of radiation field intensities, with [U^/^^, Vmaxl ranges of 
[3.8-4.8] for Region 6, [3.0-4.0] for Region 23. As observed previ- 
ously, the PAH fraction is very low in the H ii regions compared to 
the diff'use ISM around those regions. 

Finally, we also estimate the Hi mass in each aperture. Val- 
ues are tabulated in Table [2] As mentioned previously, the value 
of the Galactic G/D is 157. The LMC having a lower metallicity, 
the G/D of the LMC should be higher than this value. We observe 
a large spread in the H i/D values from one region to another, and 
low ratios in H ii regions, which is not surprising. Indeed, the whole 
region being a very large molecular complex, the molecular gas is 
thus expected to be a significant fraction of the total gas mass in 
those regions as well. As previously mentioned, further investiga- 
tions including CO observations for the whole complex will be in- 
vestigated in a future paper. 



5.3 SPIRE/LABOCA colour-colour diagram 

Figure [TT] shows Ssso/Sgvo versus S250/S350. We remove the radio 
contamination from the LAB OCA flux densities used in this plot. 
The diff'erent lines are colours obtained modelling single temper- 
ature modified blackbodies with an emissivity index of 2.5, 2, 1.5 
and 1 and generated with temperatures ranging from 5 (left) to 70K 
(right). The bright star-forming regions have higher S250/S350, con- 
sistent with regions with higher temperature dust. They also show 
lower S350/S870 ratios. N158 (region 3), N159 (region 1) and N160 
(region 2) thus present "flatter" submm slope than regions of the 
diffiise ISM (i.e. higher 870 yt/m flux density for a fixed 350 fim 
flux density). Region 5 does not appear in the diagram because of 
its very high 350-to-870 colour. Figure [9] (bottom panel) suggests 
that the 870 fim flux density we derive could be underestimated 
compared to what is expected from the SPIRE bands. In Fig. [TT] 
we indicate by a red cross where the region would fall in the plot if 
the modelled 870 fim value (derived from our fit) was used instead 
of the flux density derived directly from the 870 jdm map. 
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Figure 11. SPIRE/LABOCA colour-colour diagram. The 350-to-870 /urn 
flux density ratio are plotted as a function of the 250-to-350 //m flux density 
ratio. The brightest star-forming regions of the complex (1 to 6) appear 
in red. The individual subregions of N 159 (la, lb and Ic) are in orange. 
Regions 7 to 24 are in green. The lines show the colours obtained with 
single temperature modified blackbodies with an emissivity index of 2.5, 
2, 1.5 and 1. Radio contamination has been removed from the 870 //m flux 
densities. The red cross indicates where region 5 would fall if the modelled 
870 jum value was used instead of the flux density we estimate. 



5.4 About submm excess 

Ground-based data and now Herschel observations of low- 
metallicity environments at submm wavelengths have helped us 
to better probe the cold phases of dust and investigate variations 
in their properties. These observations have led to the detection 
of a "submm excess" ( (Galliano et al.|[200 F Marle au et al.||2006l 
[Galametz et al.||2009] ), namely submm emission higher than that 
predicted by FIR data and standard dust models ( [Praine & Li|2007| 
for instance). Several interpretations have been analysed to explain 
this excess. It could be linked with a diff'erent dust component: 
very cold dust (e.g. Galliano et al. 2003 2005, Galametz et al.] 
[2010| ), emission from "spinning dust" ( Y sard & Verstraete,2010j) 
or magnetic dipole radiation from magnetic nano-particles ( [Draine[ 
[& Hensley[2012| ). The excess could also be linked with a modifica- 
tion of the dust emissivity properties at colder temperatures ( [Meny[ 
[et al.|2007|[Paradis et al.|20T0| . 

Several studies have been dedicated to the investigation of 
submm and millimeter excess in the LMC. [BotetaTl ( [2010b| ) found 
excess millimeter emission while modelling the integrated SED of 
the LMC using FIR and radio data. This excess was explained by 
possible emission from "spinning dust". However, this hypothe- 
sis could require extreme excitation conditions for the small dust 
grains expected to be the driver of such emission. Recent Planck 
observations of the LMC were analysed in [Planck Collaboration| 
( 2011 ) and also lead to the detection of a millimeter excess long- 
ward 500 yum. They found that CMB fluctuations would be suffi- 
cient to explain the excess. Using Herschel observations of a large 
strip across the LMC, [Gordon et alT] ( [2010| ) showed, on a more re- 
solved scale, that the Herschel 500 jim emission was ~10% higher 
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than expected from extrapolating models fitted at shorter wave- 
lengths (and using a modified blackbody with an emissivity index 
of /?=1.5) and found a weak anti-correlation of the excess with the 
MIPS 24 yL/m emission or the total gas surface density. Using the 
resolved SED modelling technique we use in this paper (and AC 
to model carbon dust), [Galliano et al.| ( |2011| ) confirm the detec- 
tion of a 500 yL/m extended emission excess in the same strip across 
the LMC, with an average relative amplitude of 15 % (and vari- 
ations up to 40 %). They found that this resolved submm excess 
anti-correlates well with the dust mass surface density and proba- 
bly not linked with very cold dust or CMB fluctuations. 

The SEDs presented in Fig.|9]and Fig.[TO|do not show excess 
emission long-ward 500 jum. Our model seems sufficient to prop- 
erly fit the observations in the strongly star-forming regions as well 
as in the more diffuse regions we modelled, without requiring any 
extra component to explain the submm emission. We nevertheless 
note that unlike Gord on et al.| ( |2010l ) or [Galliano et a l. ( 2011 ), we 
do not study regions with very low surface densities in this anal- 
ysis. If the environment where submm excess is hiding is in the 
low surface brightness regions, we will not be able to detect it in 
our star-forming complex. Low surface brightness structures are 
also the regions where our data reduction of LAB OCA data is the 
most likely to have removed signal, which could explain the non- 
detection. 



6 CONCLUSIONS 

In this paper, we combine Spitzer Space Telescope, Herschel Space 
Observatory and LAB OCA mid-infrared to submm data (from 3.6 
to 870 jim) of the resolved massive star-forming region N158- 
N159-N160 located in the Large Magellanic Cloud. 

• Using a resolved 2MBB fitting technique and a Pc fixed to 1.5, 
we find an average temperature of ~27K across the region, with a 
maximum in N160. Using data from 3.6 to 500 fim, a physical SED 
modelling and the hypothesis of amorphous carbon to model car- 
bon dust, we also derive maps of the dust temperature distribution, 
the average starlight intensity, the PAH fraction, the star formation 
rate and the dust mass surface densities. The PAH fraction strongly 
decreases in our Hn regions. This decrease coincides with peaks 
in the mean radiation field intensity map, which is consistent with 
the well-known destruction of PAHs by strong radiation fields. The 
absence of PAHs could also be linked with the young age of the 
clusters, namely that the regions could be too young to have formed 
any PAH. 

• We analyse how submm colours vary across the star-forming 
complex and study their dependency with the surface brightnesses 
or the radiation field intensity <U>. We find that the submm ra- 
tios correlate strongly with <U> and that while the 100/160 and 
160/250 colours correlate more strongly with the 24 jjm bright- 
ness, the 250/350 and 350/500 colours equally correlate with the 
3.6 and the 24 fim surface brightnesses. Possible variation of the 
emissivity with temperature have been investigated using SPIRE 
colour-colour diagrams. No clear trend is observed. 

• The dust surface densities follow the FIR distribution and a 
non-negligible amount of dust is detected in N159S. We find a total 
dust mass of 2.1 x 10^ M© in the resolved elements we model (and 
2.8 times more dust if graphite is used to model carbon dust). 

• Considering that gas and dust are closely mixed in the N159 
region, we use our dust mass map to compare the distribution of 
dust with those of the tracers of the atomic gas in the complex to 
investigate variations in the atomic +molecular gas-to-dust ratio in 



N159. We use a CO(3-2) map of N159 to quantify the molecular 
reservoir in that region. A mean value of ~356 is derived for the 
N159 complex when using the Xco conversion factor of Fukui etaL| 
(2009) (=7 X 10^0 H2 cm-2 (K km s'^-^. If we consider a constant 
G/D in the complex, we can apply the minimization technique of 
[Sandstrom et al.|(2012b| ) to inversely derive the Xco conversion 
factor of the complex that minimises the D/G scatter in N159. The 
D/G scatter is minimized when the conversion factor Xco is equal 
to 8.6 Mo pc-2 (K km ^'^y^ (=5.4x10^0 H2 cm-^ (K km s'^"^). 

• We finally model the SEDs of 24 selected regions, now includ- 
ing the 870 fim data in the fitting, and describe the variations in the 
dust thermal emission (thus in the SED shape) across the complex. 
We show that our 'AC" model is sufficient to explain the submm 
emission we observe at 870 fim. 
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Table 2. Cold dust temperatures, dust masses, PAH fraction, radiation field intensity parameters and H i-to-dust mass ratios obtained for our selected regions 
using the "AC model" 



Region Mdust ^PAH a ^min ^max Hi/D 

(K) (Mo) 



1 (N159) 


34.7 


3.08 X 10^ (8.67 X lO^f 


0.37 


2.38 


6.9 


1.2 X 10^ 


165 


la (N159E) 


27.7 


8.01 X 10^ (3.40 X 10^) 


0.27 


1.94 


2.2 


3.5 X 10^ 


55 


lb 


34.7 


1.71 X 10^ (6.03 X 10^) 


0.23 


2.26 


8.2 


1.0 X 10^ 


79 


Ic (N159W) 


31.0 


7.50 X 10^ (2.92 X 10^) 


0.41 


2.32 


5.1 


3.9 X 10"^ 


67 


2 (N160) 


33.7 


1.15 X 10^ (4.32 X 10^) 


0.31 


2.27 


7.1 


1.1 X 10^ 


102 


3 (N158) 


33.6 


7.21 X 10^ (2.17 X 10^) 


0.32 


2.31 


6.5 


9.4 X 10"^ 


130 


4 


26.8 


5.29 X 10^ (1.44 X 10^) 


0.63 


2.50 


2.8 


2.2 X 10^ 


199 


5 


34.7 


5.14 X 10^ (1.36 X 10^) 


0.39 


2.39 


4.6 


1.0 X 10^ 


195 


6 (N159S) 


24.8 


6.44 X 10^ (1.92 X 10^) 


0.83 


2.42 


3.8 


4.8 


212 


7 


28.7 


3.01 X 10^ (8.65 X 10^) 


0.60 


2.50 


3.4 


3.1 X 10^ 


99 


8 


27.2 


2.65 xlO^ (9.18 xl02) 


0.54 


2.24 


2.7 


3.9 X 10^ 


131 


9 


26.7 


1.66 X 10^ (4.91 X 10^) 


0.53 


2.32 


2.8 


2.1 X 10^ 


111 


10 


29.7 


1.34 X 10^ (4.04 X 10^) 


0.44 


2.36 


4.4 


1.3 X 10^ 


126 


11 


28.5 


1.42 X 10^ (3.87 X 10^) 


0.60 


2.43 


3.2 


1.0 X 10^ 


151 


12 


26.8 


1.08 X 10^ (2.47 X 10^) 


0.86 


2.50 


2.3 


5.5 X 10^ 


164 


13 


28.2 


1.29 X 10^ (3.10 X 10^) 


0.48 


2.50 


2.9 


4.8 X 10^ 


205 


14 


30.3 


1.45 xl02 (3.66 xlO^) 


0.47 


2.50 


4.3 


4.2 X 10"^ 


196 


15 


26.8 


7.34 X 10^ (2.16 X 10^) 


0.54 


2.18 


2.7 


9.0 X 10^ 


374 


16 


24.9 


1.34 xl02 (3.57 xl02) 


0.67 


2.50 


2.0 


2.0 X 10"^ 


138 


17 


27.5 


1.84 X 10^ (4.60 X 10^) 


0.62 


2.50 


2.8 


4.1 X 10^ 


145 


18 


27.2 


9.20 X 10^ (2.25 X 10^) 


0.66 


2.50 


2.4 


8.2 X 10^ 


219 


19 


24.2 


1.09 X 10^ (3.30 X 10^) 


0.48 


1.75 


1.0 


7.3 X 10^ 


202 


20 


24.5 


1.64 X 10^ (5.37 X 10^) 


0.78 


1.78 


1.1 


5.5 X 10^ 


176 


21 


25.6 


8.79 X 10^ (3.23 x 10^) 


0.74 


1.00 


1.0 


2.1 X 10^ 


327 


22 


24.8 


9.97 X 10^ (2.47 x 10^) 


0.67 


2.50 


1.9 


3.5 X 10^ 


275 


23 


24.1 


1.07 X 10^ (2.60 X 10^) 


0.95 


2.50 


3.0 


4.0 


289 


24 


23.9 


7.76 X 10^ (1.92 X 10^) 


0.99 


2.50 


1.3 


4.6 X 10^ 


182 



^ We indicate the dust masses obtained using the graphite model in parenthesis for comparison. 



